Universal Aspects of Coulomb Frustrated Phase Separation 
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We study the consequences of Coulomb interactions on a system undergoing a putative first order 
phase transition. In two dimensions (2D), near the critical density, the system is universally unstable 
to the formation of new intermediate phases, which we call "electronic microemulsion phases," which 
consist of an intermediate scale mixture of regions of the two competing phases. A correlary is that 
there can be no direct transition as a function of density from a 2D Wigner crystal to a uniform 
electron liquid. In 3D, if the strength of the Coulomb interactions exceeds a critical value, no phase 
separation occurs, while for weaker Coulomb strength, electronic microemulsions are inevitable. 
This tendency is considerably more pronounced in anisotropic (quasi 2D or quasi ID) systems, 
where a devil's staircase of transitions is possible. 



We consider the effect of long range Coulomb interac- 
tions on a system undergoing a first order phase tran- 
sition between two compressible states. In the absence 
of Coulomb interactions, a first order transition implies 
an interval of mean density in which the equilibrium 
state is macroscopically phase separated into regions of 
higher and lower than average density. A Coulomb in- 
teraction precludes macroscopic phase separation; con- 
sequently, the system can either undergo a direct, first 
order phase transition between the two competing uni- 
form states at the critical density, or can form interme- 
diate phases, which we refer to as "electronic microemu- 
sions," which can be thought of as a mesoscale mixture 
of the two competing phases. It has been suggested that 
such states resulting from Coulomb frustrated electronic 
phase separation are relevant to the physics of various 
highly correlated materials, including the cuprate high 
temperature superconductors and the collosal magneto- 
resistance manganatesQ], 0, as well as toy, Q the pure 
two dimensional electron liquid at large r s (low density) . 

In the present paper, we consider several universal (i.e. 
independent of microscopic details) aspects of Coulomb 
frustrated phase separation. In two dimensions (2D), we 
show that even in the absence of any quenched disorder, a 
direct first order phase transition as a function of density 
between two distinct thermodynamic phases is forbidden. 
Rather, in the neighborhood of the putative critical den- 
sity, there is an instability of the uniform state to some 
form of mesoscale phase-coexistence, leading to new in- 
termediate phases. In particular, this theorem applies 
to the 2D Wigner crystal to electron fluid transition. 

In 3D, a direct first order transition (with no region of 
phase coexistence) occurs if the strength of the Coulomb 
interaction, Q, exceeds a critical^ strength, Q > Q c , 
while if Q < Qc, one or more microemulsion phases oc- 
cur. Under many circumstances, at least at mean-field 
level, the transition from homogeneous to microemul- 
sion phases as a function of decreasing Q is a Lifshitz 
transition, in which the characteristic period of the mix- 
ture diverges as Q — > Q c from below, and there can be 
a Devil's staircase consisting of a complex mixture of 



commensurate and incommensurate density waves. For 
Q "C Qc, the characteristic size of the regions of each 
phase is small, and more detailed, microscopic consid- 
erations become important. In anisotropic (quasi 2D or 
quasi ID) systems, Q c becomes increasingly large with in- 
creasing anisotropy, making some form of Coulomb frus- 
trated mesoscopic phase separation nearly inevitable. 

An explicit model: The analysis described in the 
present paper is quite general, but it is nevertheless more 
convenient, pedagogically, to present it in the context of a 
simple model described by the following classical Hamil- 
tonian density 



n = E t- [M 2 + u(<p) -fx(p- Pc )<t> + 



a = l 



-[p-p]V[ P -p], 



(1) 



where U is a potential, which to be concrete we take 
to be U = \[<f) 2 — l] 2 /4, p is the charge density, p is the 
(non-dynamical) uniform background charge density, and 
we have expressed the Coulomb interaction, in an oper- 
ator form, such that Vf = J d r'V(f — r')f(f ) with 
V(r) — Q/r. In this model, the two uniform phases have 
<\> = ±1 + p(p — p c )/2X + . . . which we label as 4> w ±1, 
assuming A to be large. The term in TL proportional to 
p expresses the fact that the (t w 1 phase is favored at 
higher densities, and the <f> « — 1 phase is favored at low. 
For an isotropic 3D system, while we will call a 

system quasi-2D if there exists one direction in which K a 
is very small, and quasi ID if there are two such direc- 
tions. (We will label spatial directions in such a way that 
Ki > ft 2 — K 3-) We- use a uniform continuum notation, 
but where discreteness is important (such as in a layered 
system, in which position in the perpendicular direction 
is labeled by a layer index), it is to be understood that 
d a 4> = [4'{ n + 1) — 4>{n)]/a a where a a is the lattice con- 
stant. The terms . . . represent higher derivative terms, 
and terms proportional to higher powers of p — p. 

The significance of n a and A is that they determine 
the local (i.e. the Q — > limit) structure of an in- 
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terface between the two phases running perpendicular 
to the a direction: Specifically, the surface tension is 
<r a = (7t/2)a/Ak q ,/2 and the width of the interface is 
a a = y K a /2X. 

Because the Hamiltonian density is a quadratic func- 
tion of p, density fluctuations can be integrated out to 
yield an effective Hamitlonian for <f>, alone. The result 
can be expressed formally as 

D 

ueff = ^ >^ M2 + u{(f>) _ _ pc)(f) + 



1 

- ^2 d a <l>Vd a (j), 



(2) 



a=l 



where the renormalized interaction, V , is defined in terms 
of the inverse of V according to — V 2 V = p 2 V^ x , or 
taking the Fourier transform: Vg = /i 2 [fc 2 Vg] _1 . 

The 2D case: In 2D, the Fourier transform of V is 
Vg = p 2 [TrQ\k\]~ 1 . Transforming back to real space, the 
final term in Eq. yjis 



E. 



Coulomb 



4ir 2 Q 



(3) 



In the limit in which the interfaces between phases are 
narrow compared to the size of domains, one can neglect 
the gradient of <f>(r) in the bulk of the phases, so the 
expression for the total energy can be simplified to read 



E = n[p-Pc}[A+ -AJ\ 



dla{0) 



dl ■ dl' 



K 2 QJ \l-l'\ 



(4) 



where the integrals run along the interface in the direc- 
tion such that the (f> « +1 phase is to the left, A± are the 
areas of the two phases, a(9) is the microscopic contri- 
bution to the surface tension, which can depend on the 
direction of the interface, and any short-distance singu- 
larities in the final integral are cutoff at the scale of the 
interface width, a. 

Eqs. and 01 are the principle results of our paper; the 
classical low temperature thermodynamics is obtained 
by integrating over all domain patterns with Boltzman 
weight determined by this effective Hamiltonian. The 
key feature of this sum is that the third term in Eq. 0] 
is negative and logarithmically divergent whenever the a 
domain is large. Thus, close enough to the point of the 
phase transition, micro-emulsion phases which are phase 
separated on a mesoscopic scale always have lower en- 
ergy than the uniform phases. Q It has previously been 
shown0 on quite general grounds that if there is a direct 
transition from a crystal to a uniform fluid, it must be 
first order. Combined with the present result, this implies 
that there can be no direct transition. Rather, there must 
exist one or more intermediate phase, and a sequence two 
or more continuous transitions [lfj. Moreover, in the case 
of the the 2DEG, the first of these transitions must occur 



at strictly higher density (smaller r s ) than the putative 
Wigner crystal to uniform liquid transition. |16| 

To find the mean-field phase diagram one should mini- 
mize Eq. Uwith respect to the geometry of the interfaces. 
Remarkably, the same expression for the interfacial en- 
ergies (the third term in Eq. 0} arises 0, 0, 0, 0] in 
a 2D problem where phase separation is frustrated by 
long-range dipolar interactions (although in this case the 
first, bulk contribution to the energy is different). Thus 
what is known about that problem can be easily applied 
to the present case. 

On the mean field level the shape of the micro-emulsion 
depends of a degree of anysotropy of the surface tension. 
The simplest situation arises if a is highly anisotropic - in 
this case, the domain walls always lie preferentially along 
the easy direction and thus, at least at mean-field level, 
the lowest energy microemulsion phases are striped Q . 

The energy per unit area, relative to the uniform phase, 
of an alternating array of stripes of <f> « +1 of width W 
and <f> w — 1 of width L — W is computable from Eq. 0] 

e = 2L~ l {p Sp W + cr[l -7log(W/o) -jf(W/L)]}(5) 

where Sp — p — p c , 7 = p 2 /(ir 2 aQ), and f[x) = 
log[sm(7rx)/7ra;]. For p — p Cl this energy is minimized 
by alternating stripes of equal width, L = 2W, with 
W = (2/tt)W and 



W = aexp[l + I/7] 



(6) 



is a characteristic emergent length scale. For p > p c , 
the ratio of Wj L decreases monotonically, until as Sp — > 
Ap = (ctj/p^Wq, the period of the stripes diverges as 



W i 



/ 3C(2)Ap 
Ap-Sp 



(7) 



and the width of the stripes, W — ► Wo, where ((2) = 
7r 2 /6 is the zeta function. 

Thus, as a function of decreasing p, there is a "Lif- 
shitz" transition jl.8] from the uniform phase to the stripe 
phase at p = p c + Ap, then a continuous evolution of the 
stripe phase with the identity of the minority and major- 
ity phases interchanging at p — p Cl and finally a second 
Lifshitz transition at p = p c — Ap. The resulting mean 
feild phase diagram is qualitatively identical to the phase 
diagram in the case with dipolar interactions - see Fig. 
2c of Ref. However, in the dipole case, the density 
differences between the two phases is roughly determined 
by a Maxwell construction, and the characteristic stripe 
width is determined by the size of the dipole, d, which is 
an independent physical parameter, which in some phys- 
ical realizations can be large compared to the spacing 
between particles. In the present case, the density con- 
trast and the characteristic width of the stripes, Wo, is 
determined by microscopic physics. 

When the surface tension is more nearly isotropic, the 
mean-field phase diagram is more complicated. Although 
for p near p c , stripe phases continue to have the lowest 
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energy, near the transition to the uniform phase, bub- 
ble phases, consisting of isolated bubbles of the minority 
phase, have lower energy than the stripe phase. The Lif- 
shitz points at p = p c ± Ap now mark the transition 
between a uniform and a long-period bubble phase 0, Q . 
In the case of the triangular Wigner crystal-Fermi liquid 
transition the anisotropy of the surface tension is such as 
to produce stable bubble phases. Although the optimal 
bubble is facetted, one can establish the stability of the 
bubble phase by using as a variational ansatz hexaganol 
bubbles. On approach to the Lifshitz point, the size of 
the bubbles is of order Wo while the distance between 
them diverges as 

L B ~W [Ap/(Ap-6p)]-k (8) 

This is not quite the end of the story; since the transition 
between the bubble and stripe phases is first order, a sec- 
ond level of frustrated phase separation occurs producing 
a new set of intermediate phases. These states consist of 
alternating stripes of stripes and stripes of bubbles. The 
resulting phase diagram is qualitative identical to that 
shown in Fig. 2d of Ref. [4( for the 2D dipolar case. 
Again, the difference is that in the dipolar case the size 
of the bubbles is proportional to the size of the dipole, 
and can be parametrically large. 

Because of the importance of the 2D case, it is worth 
deriving Eq. 0|in a second way. For simplicity, consider 
the situation in which the background charge density is 
tuned to the critical value at which the putative first 
order transition occurs, p — p c . We compute the dif- 
ference in energy between the uniform density state and 
that in which there is an interface such that (other than 
in a narrow interface region of width 2a), there is phase 
4> w —1 for x < —a and phase <j) ~ +1 for x > a. So as 
to minimize the Coulomb cost of this interface, while at 
the same time gaining maximal energy from phase sepa- 
ration, we allow the charge density to assume the profile 
p{r) = p + (j>(x)A(a\x\)^ 1 . (This expression is identical 
to that arising in the problem of a contact potential of 
two metals with different work functions 0.) With this 
profile, we easily see that the Coulomb cost of the in- 
terface is E c — A 2 QL\og(W/a), where L is the length 
of interface and W is the transverse width of the two 
regions separated by it. Similarly, the energy of phase 
separation gained within each region by the density re- 
distribution is Eps = — 2pA(L/a) \og(W/a). Because 
the major contributions to these logarithms comes from 
regions far from the interface, where the deviations from 
uniform density are small and slowly varying, the result 
is insensitive to the inclussion of any further terms in the 
energy functional. The energy of this interface is mini- 
mized by A = p/aQ, at which point 

E in terface = [oo - (p 2 /aQ) log(W/a)]L (9) 

which is manifestly negative for large enough W. 

Taking, as an example, a striped state, where W is the 
width of the stripes, the logarithm in Eq. [5] can be seen 



to be precisely the same one that comes from integrating 
the expression in Eq. 21 over the interface, Eq. [5] Eq. 
E3 by itself, constitutes a proof that a first order phase 
transition is forbidden - it implies an absolute instabil- 
ity of any uniform (or macroscopically phase separated) 
state at the point of the putative transition. Although 
the analysis that leads to this conclussion is classical, be- 
cause it is a long-wave-length instability, it is unaffected 
by quantum fluctuations, and so applies at zero temper- 
ature, T = 0, as well as non-zero T. 

The characteristic size of bubbles and stripes depends 
exponentially on the ratio of microscopic energies, 7: If 
7 is roughly 1 or more, then the size of the domains 
is microscopic. Moreover, in this case, there is no rea- 
son to expect quantum or thermal fluctuations to be 
small. Consequently, although the instability of the uni- 
form phases is a robust long-distance phenomenon, the 
mean-field phases we have found may or may not survive 
these fluctuations. However, because of the exponential 
dependence, it can happen that 7 is small compared to 
1, in which case Wq can be very large compared to a. In 
this case, the mean- field analysis presented here should 
be reliable, and Wq is a non-trivial emergent length. The 
significance of the classical and quantum fluctuations has 
been discussed in 0, • 

The 3D case: In 3D, it follows simply from the fact 
that the Fourier transform of the Coulomb interaction is 
V(k) = 4nQ/k 2 , H eff is local 

neff = £ K - M 2 /4^Q] [d ^ ? + m _ m + 

(10) 

Thus, if the Coulomb interactions are sufficiently strong 
that Q > A c = /i 2 /47TK 3 , then as a function of p, the 
system undergoes a direct first order transition from a 
uniform <f> ~ 1 phase for p > p c , to a uniform <f> w — 1 
phase for p < p c . In contrast to 2D, there is no absolute 
proscription against first order transitions in 3D. 

However, for weaker Coulomb interactions, for p in the 
neighborhood of p c , there is an intermediate modulated 
phase whose precise structure is determined by the higher 
derivative terms which are not explicitly exhibited. (In 
part, this case has been previously considered in Q.) 
For illustration, consider the phase diagram at the crit- 
ical density, p — p c . At the critical value, Q = Q Cl the 
coefficient of the leading stiffness term changes sign, the 
uniform state must be unstable, and higher order elas- 
tic constants (more derivatives) must be included in the 
effective Hamiltonian. Most simply, we can include the 
next order terms in Eq. ^3 ■ ■ ■ = J2 a K a(^a^) 2 + ■ ■ • • So 
long as these higher order elastic constants remain posi- 
tive, the transition to the modulated phase as a function 
of Q is the classic Lifshitz transition. Specifically, as Q 
approaches the critical value Q c from below, the period of 
the modulated phase diverges as L ~ a^J Q c / (Q c — Q). 
Alternatively, if k' 3 < 0, this Lifshitz transition is pre- 
empted by a first order transition (at a larger critical Q) 
to a modulated phase with a short period. 
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Layered systems: A natural realization of a quasi 2D 
system is in a layered material. This is still a 3D system, 
and so subject to the analysis of the previous paragraphs, 
but now there can be significant effects of the discrete, 
lattice structure. In particular, so long as the spacing 
between layers, 03, is not too large, this discreteness in- 
troduces the only important modification of the above 
discussion. Where the continuum theory predicts a dis- 
continuous transition to the modulated phases, lattice 
discreteness generally locks the period of the resulting 
modulation to a specific, low-order commensurate struc- 
ture. Alternatively, near the Lifshitz transition, the com- 
petition between the period favored by the elastic con- 
stants and the underlying lattice typically results in a 
Devil's staircase of modulated phases, and moreover can 
lead to intrinsically glassy behavior associated with the 
pinning of domain walls. 

The easiest way to get a flavor of the resulting pos- 
sibilities is by simply discretizing the spatial derivatives 
in Eq. 1101 in the case in which n\ — K2 > K3, and 
n' a > 0, the resulting model is the soft-spin version 
of the much studied 13] antiferromagnetic next-nearest- 
neighbor Ising (ANNNI) model, in which the antiferro- 
magnetic coupling between second neighbor planes is pro- 
portional to K3, the couplings between nearest-neighbor 
planes has a strength (and sign) which varies as a func- 
tion of K3 — [i 2 /4irQ, and the in plane nearest-neighbor 
ferromagnetic couplings are proportional to K\. This 
model has a remarkably subtle and beautiful phase dia- 
gram consisting of uniform (ferromagnetic) phases, short 



period commensurate phases, of which the most promi- 
nent are the period 2 (alternating up and down planes) 
and period 4 (2 up 2 down) phases, and then a "Devil's 
flower" consisting of high order commensurate, and in- 
commensurate phases. Presumably, the more subtle as- 
pects of this phase diagram are lost when thermal or 
quantum fluctuations are included, but the tendency for 
Coulomb frustrated phase separation to produce patterns 
of alternating high and low density layers is robust. 

In the limit of larger 03, the phase diagram is still 
more complex. So long as 03 is less than the character- 
istic size, Wq, which characterizes the domain size in an 
isolated layer, interlayer phase separation is energetically 
preferred over intra-layer phase separation. However, for 
a 3 3> Wo , 2D patterns of phase separation dominate the 
physics at shorter length scales, while intralayer consider- 
ations become important only at much longer distances. 

Quenched disorder produces its own form of locally 
pinned mixtures of the two competing phases. Distin- 
guishing the two effects in experiment requires studying 
the properties of progressively cleaner systems. With this 
caveat, we believe the present results are significant for 
a host of phe nomena in layered and quasi-lD crystalline 
materials [T^. as well as for0,0| the 2DEG in high mo- 
bility semiconductor devices at large r s . 
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